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ABSTRACT 

We report a new Fabry-Perot search for Ha emission from the intergalactic 
cloud HI 1225+01 in an attempt to measure the low redshift ionizing background 
radiation. We set a new 2a upper limit on Ha emission of 8 mR (5 x 10~^° 
ergs cm~^ s~^ arcsec"^). Conversion of this limit to limits on the strength of 
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the ionizing background requires knowledge of the ratio of the projected to total 
surface area of this cloud, which is uncertain. We discuss the plausible range of 
this ratio, and within this range find that the strength of the ionizing backround is 
in the lower range of, but consistent with, previous observational and theoretical 
estimates. 

Subject headings: Cosmology: diffuse radiation. Galaxies: intergalactic medium 

1. INTRODUCTION 

The value of the local ionizing fiux is of interest for several reasons: 1) It plays an 
important role in setting the level of ionization of the low density hydrogen gas which gives 
rise to the Lja forest. Recent simulations have been remarkably successful in reproducing 
the main features of the Lya forest from high redshifts to zero redshift (Dave & Tripp 2001). 
As discussed in detail in §4, for a given set of cosmological parameters, these simulations 
are now able to make predictions about the value of the local H I photoionization rate, 
and it would be valuable to determine empirically the value of this rate to compare it 
with the simulations. With a sufficiently accurate determination of its value and further 
forthcoming data on the low redshift Lya forest, together with further improvements in the 
simulations it may be possible to set some constraints on these cosmological parameters, 
especially the baryon number. 2) A limit on the photoionization rate in locations which 
are free of local ionizing sources (i.e. star-forming regions or ionization due to high energy 
particles in the neighborhood of galaxies) can then be used to assess the relative importance 
of such local sources near a galaxy compared to the non-local, global value. In particular, 
there is currently some controversy concerning the ionization source for some high velocity 
clouds in our own galaxy (Weiner, Vogel, & Williams 2000; Bland-Hawthorn & Putnam 
2001). 3) Estimates of the photoionization rate due to the cumulative effect of active galactic 
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nuclei can now be carried out with fair confidence (Shull et al. 1999, hereafter SRGPF); but 
there remain significant uncertainties associated primarily with the luminosity functions of 
these objects and with the evaluation of the opacity associated with intervening neutral 
hydrogen. 4) Similar remarks apply to the contribution from escaping radiation shortward 
of the Lyman limit from hot stars in galaxies. However, this estimate is also dependent upon 
the fraction of the Lyman continuum radiation assumed to escape from the galaxies, and 
direct observational evidence about the magnitude of this escape fraction is still somewhat 
scanty. 

For these reasons, we have previously set limits on the value of the non-local ionization 
rate at zero redshift by setting limits on the Ha flux from an isolated, large, H I cloud 
which is optically thick to the Lyman continuum (Vogel et al. 1995, hereafter VWRH). 
The essence of the method is very simple. Illumination by ionizing radiation of an H I 
cloud that is optically thick to the Lyman continuum results in Ha emission that is well 
described by Case B recombination theory. Accordingly, approximately one Ha photon 
escapes for every two incident Lyman continuum photons; the exact value depends only 
mildly on parameters that are reasonably well known. The one significant uncertainty is the 
ratio of the observed projected surface area to that illuminated by the ionizing radiation; 
however, the uncertainty is perhaps a factor of two, and can be better determined for objects 
where the geometry is known, such as galactic disks. The expected Ha radiation is extremely 
faint — based on predictions for the low redshift metagalactic ionizing flux (e.g., SRGPF), 
the Ha surface brightness is expected to be only a few milliRayleighs (for Ha, 1 mR = 
5.67 X 10^^^ ergs cm^^ s^^ arcsec^^). By comparison, the observed Ha emission from objects 
ionized by other sources is significantly brighter: high velocity clouds known to be within 
a few kpc are typically 60-200 mR (Tufte, Reynolds, & Haffner 1998); most compact high 
velocity clouds lie between 30 to 300 mR (Weiner, Vogel, & Williams 2000); the Magellanic 
Stream ranges up to more than 1000 mR (Weiner & Williams 1996); a brightness of 90 mR 
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was measured in the outer disk of NGC 253 (Bland-Hawthorn, Freeman, & Quinn 1997); and 
the diffuse extraplanar ionized gas imaged in edge-on spirals is typically brighter than 1000 
mR (Rand, Kulkarni, & Hester 1990; Veilleux, Cecil, & Bland-Hawthorn 1995). As a result 
of its expected faintness, detection of Ha emission due to ionization by the metagalactic 
ionization requires observation of extended clouds using specialized techniques. 

There exist other approaches to an empirical value for the global zero redshift ioniza- 
tion rate. However, several of these methods suffer from significant limitations. At higher 
redshifts, the proximity effect (Bajtlik, Duncan, & Ostriker 1988; Bechtold 1994) has been 
used to estimate the value of the metagalactic ionizing flux; however at low redshift the 
sparseness of higher resolution data and the relative sparseness of the Lya forest itself makes 
this method currently less useful. In addition, there are some issues about the interpretation 
of the proximity effect (e.g., the clustering properties of matter in the vicinity of QSOs; see 
Scott et al. (2000) for further discussion). 

A potentially very powerful technique involves observing the rapid transition from a 
mostly neutral hydrogen disk to a mostly ionized disk as one moves outward along the 
radius of spiral galaxies. From the value of this transition column density, one can infer 
the flux of the background ionizing radiation given a model for the vertical distribution of 
gas. The basic idea for this technique goes back to papers by Sunyaev (1969) and Felten & 
Bergeron (1969). The historical development of this idea is summarized by Maloney (1993) 
who applied it to realistic galaxy models. Maloney gives a careful treatment of the radiative 
transfer equation and discusses a wide range of galaxy parameters. His results show that 
the derived ionizing flux is relatively insensitive to the dynamical models adopted for the 
galaxy. The application to NGC3198, which has the most sensitive HI observations in this 
context, and whose geometry and dynamics are not very uncertain, yielded a transition 
column density of about 5 x 10^^ cm"^ and from this Maloney deduced that the flux of 
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ionizing photons is in the range 5 x 10^ to 5 x 10^ photons cm~^ sec~^. Shortly after this 
work Dove & ShuU (1994) pubhshed a similar analysis and confirmed these results. They also 
showed that the results are not very sensitive to the assumed spectral index of the ionizing 
background. 

Application of this method to other objects has also been carried out. In particular, as 
noted above, measurements of Ha emission beyond the HI truncation radius were obtained 
by Bland-Hawthorn, Freeman, & Quinn (1997) in NGC253. In this case, however, these 
authors noted that the inferred surface brightness was higher than could likely arise only 
from the metagalactic radiation field. In addition, the ratio of [Nil] A6548 to Ha emission 
was unusually high. These authors favored a model in which a warp in the outer disk allowed 
ionizating flux from the inner young star population to dominate the ionization. In NGC253 
the ratio of the HI disk radius to the optical disk radius is considerably smaller than in 
NGC3198. Nevertheless, Dove & Shull (1994) do not reject the possibility that hidden local 
sources of ionization may contribute significantly to the ionizing flux even in NGC3198. 

We note that there appear to be instances where no abrupt truncation of H I occurs to 
column densities which are smaller than ~ 5 x 10^^ cm^^. For example, Walsh, Staveley- 
Smith & Oosterloo (1997) present HI images of NGC 289, a giant gas-rich low-surface- 
brightness spiral, and show HI contours down to a column density of 1 x 10^^ cm~^; they note 
that no truncation is evident. In DDO 154, Carignan & Purton (1998) find no truncation, 
down to a level of 1 x 10^^ cm~^. This might be indicative of a value for the ionizing 
background somewhat lower than that derived for NGC3198, but it is also possible that the 
low observed HI column densities are a consequence of beam dilution; with better angular 
resolution, the actual column densities might be observed to be substantially higher. In 
addition, however, as emphasized by Maloney (2001, private communication), the abruptness 
of the H I cutoff depends upon the run of the total hydrogen column density which is not a 
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priori known, and the ratio of the neutral to total hydrogen column density is sensitive to the 
total hydrogen column density in the regime where the predominate ionization state switches 
from neutral to ionized (cf Figures 4 and 8 of Maloney 1993). In fact, the estimate of the 
background radiation in NGC3198 was carried out by using the observed H I column density 
to model the total hydrogen column density for assumed values of the ionizing background. 
It is the uncertainty of the true radial variation of the total column density of hydrogen at 
these large radii that appears to be the principal uncertainty in the determination of the 
background radiation using this technique. 

In any event, if the relatively large value of the metagalactic ionizing flux implied by 
the HI truncation of NGC 3198 is correct, then the Ha emission produced by this ionization 
should be detectable with achievable sensitivities in instances where the HI column density 
is high enough to ensure that most of the ionizing photons are absorbed. 

For these reasons, we continue to believe that the use of Ha surface brightness measure- 
ments as described in VWRH still provides the route to the least model-dependent estimate 
for the metagalactic ionizing flux, provided that (i) The targets are relatively isolated from 
star forming regions and (ii) that the geometry is well understood, for the reason discussed 
in section 3.3 

Consequently, when subsequent further improvements in our instrumentation were made, 
we again carried out basically the same series of observations as described in VWRH, which 
resulted in substantial further lowering of an upper limit on the zero redshift photoionization 
rate. The purpose of the present paper is to describe and discuss these observations. 

The selected target, HI 1225+01, is the same as in VWRH. This large and isolated 
cloud of H I was discovered serendipitously by Giovanelli & Haynes (1989), and subsequent 
mapping by Giovanelli, Williams, & Haynes (1991) showed it to consist of two distinct 
components — a northeast component (NE) and a southwest component (SW), with a weaker 
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bridge of H I connecting the two. The NE component has embedded in it stars and a small 
region of star formation (Turner & MacFadyen 1997, and references therein). However the 
SW component shows no trace of a stellar component (Salzer et al. 1991, and references 
therein). For this reason, we chose to observe only the SW component since the relatively 
bright Ha emission from the hot stars in the NE component would mask the weak diffuse 
signal we are attempting to detect. Subsequently, a higher resolution H I study (40") using 
the VLA was carried out by Chengalur, Giovanelli & Haynes (1995); they showed that the 
SW component appears to be a flattened inclined disk with a systematic velocity Vhei = 1260 
km s~^ and a maximum rotational velocity of ~13 km s~^. We discuss the properties of the 
SW component further in §3.3.2. The cocoon of H I in which the optical central part of the 
NE component resides makes it unlikely that ionizing radiation from the hot stars in this 
component will ionize a signiflcant portion of the SW component, and the isolation of HI 
1225+01 from other galaxies makes the SW component an attractive target for this program. 

The plan of the remainder of the paper is as follows: In §2 we describe the instru- 
mentation and observational procedures used for the results described below and in §3 we 
describe the reduction procedures. We present the results in §4 along with a discussion of 
their implications for some of the issues described above which motivated these observations. 

2. INSTRUMENTATION AND OBSERVATIONS 

2.1. Instrumentation 

2.1.1. Reducing Camera 

The observations described here were carried out at the Carnegie Observatories on Las 
Campanas, using the duPont 2.5m telescope. A reducing camera was constructed for this 
purpose. The optical design was carried out by H. Epps, with mechanical design by A. 



Schier and associates. The design consisted of a collimator with an on-axis beam size of 70 
mm and an effective focal length of 525 mm to match the f/7.5 incoming beam from the 
Cassegrain focus of the duPont telescope. A 70 mm Fabry-Perot etalon (whose properties 
are described in §2.1.2) was placed in the collimated beam. An f/2.5 camera reimages the 
telescope focal plane onto a SITe 2048 x 2048 CCD with 24 micron pixels, and a scale at the 
CCD of approximately 0.78" per pixel. The CCD was operated in a mode which yielded a 
readout noise of about 6e per pixel. However, as explained below, we were obliged to accept 
a tradeoff between adequate sampling of the signal and readout noise, and thus we binned 
the image in 4 x 4 pixels. The combined collimator-camera reducing camera is achromatic 
over a very broad range and yields an unvignetted field of about 25' in diameter. 

2.1.2. Etalon 

The University of Maryland ET-70 etalon was used for this experiment. The etalon 
and CSIOO controller are made by Queensgate. The 70 mm diameter etalon has a gap of 44 
yum, and a free spectral range (FSR) of 48.7 A at 6590.4 A, the expected center of the red- 
shifted Ha. The FWHM of the instrumental response is 1.40 A, and the effective bandwidth 
(see §3.1) is 1.64 A. The etalon was fiushed with dry nitrogen during the observations and 
calibrations. The etalon can be tilted so as to direct refiections out of the field of view. 

A four-cavity, 6-inch square, order-isolating filter was located near the telescope focal 
plane. The filter has an intrinsic FWHM of 38 A, with a transmission greater than 75% for 
6581 A < A < 6605 A. Since the etalon FSR is 48.7 A, in this wavelength range photometric 
calibration on continuum objects does not suffer from contamination from other orders. 
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2.2. Observations 

We observed HI 1225+01 the nights of 1997 April 7-10 using the reducing camera and 
Fabry-Perot etalon described in §2.1. The first three nights were photometric, while the 
atmospheric transparency varied during the course of the last night. Consequently, only 
data from the first three nights were used in the analysis described here. 

To reach the very low surface brightness required for this experiment demanded careful 
attention to observing and reduction strategy. Both shot-noise and systematic errors could 
contribute to the error budget. A key strategy for minimizing systematic errors is to switch 
as rapidly as possible between the target and sky positions. The lower limit on the switching 
time is set by the need to minimize the relative error contributed by CCD read noise. As dis- 
cussed in §2.1, the reducing camera was designed for this experiment with high throughput, 
high speed, and a large field of view. Nonetheless, due to the narrow effective bandwidth 
(1.64 A), further steps were required to make the relative contribution of CCD read- noise 
to the total error sufficiently small. First, as noted above we used 4x4 on-chip binning of 
the CCD; this yielded 3.1" pixels. Hereafter we will use "pixel" to refer to one of the 4x4 
binned pixels. Second, all exposures times were 900 seconds in length. 

To minimize errors introduced by variations in atmospheric and instrumental response, 
we adopted the following observing procedure. First, we observed blank sky at a position 
offset in RA to the west (referred to as W) of the subsequent HI 1225+01 target observation 
(referred to as GHl), with the RA offset chosen so that each pair of W and GHl exposures 
was observed through very similar azimuth and elevation ranges. In this way, we sought to 
minimize the effects of systematic errors due to azimuthal and elevation variations in the 
atmospheric, telescope, and instrument response. To eliminate linear temporal drifts, we 
repeated an exposure on HI 1225+01 (GH2) followed by an exposure on a sky position offset 
to the East (E). A W-GH1-GH2-E set of four exposures, observed in the manner described. 
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greatly reduced errors due to directional and temporal variations in the instrumental and 
atmospheric response. 

We also spatially dithered each exposure by an amount varying from 5 to 45" in order 
to minimize errors due to stars and other features. 

The monochromatic response in the focal plane is along a ring whose radius r varies 
according to the usual 

A = Ao(l + (r//f)-^/^ (1) 

where Aq is the wavelength that transmits on-axis and / is the camera focal length. As 
usual, Aq (and hence A) is adjusted by varying the etalon gap spacing, /, according to 
nA = 2 fi I cos6, where n is the order, /i is the refractive index of the nitrogen gas in the 
etalon gap, and 9 is the angle with respect to the optical axis. 

In Figure 1 we show a typical 900 second exposure toward the target position. Due to 
the large field, non-zero etalon tilt angle, and the quadratic nature of the wavelength response 
(e.g., equation 1), three to four orders are present; i.e. three to four rings are present for 
each night sky line. The on-axis order is 136, with orders 135, 134, and 133 also present 
off the optical axis. With our large, binned pixels, the quadratic nature of the wavelength 
response results in a wavelength shift comparable to the spectral resolution across a pixel for 
orders 135, 134, and 133. Since the wavelength resolution is consequently slightly degraded 
for the outer orders, we extract separate spectra for each order. 

By the well known Jacquinot relation, QR = 27r (where Q is the angular area and R 
is the resolution), the angle on the sky subtended by a complete monochromatic ring is 
independent of radius r. Clearly, we want as much as possible of the ring corresponding to 
the redshifted Ha line to be projected within the HI 1225+01 cloud. But we also want as 
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Fig. 1. — Flat-fielded 900-second exposure toward target position. The image has been 
clipped where the fiat-field response drops below 60%. Most of the rings are produced by 
OH night sky emission; the brighter lines are labeled. Note that the same lines appear in 
each order. The brightest line (OH 6577 A) is clipped in some places because the fiat-field 
response near the edge of the filter passband drops below 60%. The Galactic [N II] line also 
appears. The diameter of the field of view is ~ 25'. 

much spectral coverage on both sides of Ha as possible. Coverage at shorter wavelengths 
extends to the blocking filter cutoff at 6577 A and includes the Galactic [NH] 6584 A line and 
the atmospheric OH 6577 A line. However, to observe emission at a given A (in particular 
around 6590.4 A, the wavelength expected for the Ha centroid in HI 1225+01 SW; see §3.2), 
equation (1) shows that the maximum wavelength Aq (i.e. the wavelength observed on axis) 
increases as r is increased. In order to keep a substantial fraction of the 6590.4 A ring within 
the cloud boundary, we used a ring radius of 55-65 pixels (170-200"). This keeps roughly 
half of the ring within the cloud boundary but cuts off the red end of this spectral order 
at approximately 6593A which yields a rather short range of spectrum to the red of the 
redshifted Ha over which to define the continuum. 

The etalon and the blocking filter are normally tilted to direct refiections out of the 
field of view. To maximize the ring area on target, we chose an intermediate tilt angle for 
the etalon such that the entire ring at 6590.4 A lies within the field of view. 

As shown by equation (1), the wavelength transmitted by the etalon varies as a function 
of position. The usual method for constructing a spectrum is to obtain multiple exposures, 
stepping the etalon gap / through a full free spectral range, or at least sufficiently to fully 
sample the line profile and continuum at all positions. However, the long exposure times make 
this impractical. Instead, we know that Ha emission ionized by metagalactic radiation will 
be nearly uniform in brightness for regions that are optically thick to the Lyman continuum 
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(as described in more detail in §3.3). Consequently, we can assemble the spectrum using 
the varying spectral information provided by the different pixels in the field of view. It is 
essential that the pixels used for the Ha portion of the spectrum be limited to those pixels 
that lie towards regions of HI 1225+01 that are optically thick to the Lyman continuum. 
Pixels which transmit non-Ha wavelengths are not restricted to lie towards the cloud. 

Figure 2 shows the average sky spectra obtained for one night for orders 136 and 135, 
along with the measured transmission of the blocking filter in the tilted configuration used 
for the observations. Three bright OH atmospheric lines are evident at 6577, 6596, and 
6604 A, along with a weak OH line near 6592 A; only the 6577 and 6592 lines lie within 
the spectral coverage of order 136. Also apparent is a bump near 6583 A corresponding to 
Galactic [N II]. Continuum due to light scattered from the atmosphere, zodiacal dust, and 
Galactic sources is also present at a level of ~ 6 x 10^-*^^ ergs cm~^ s^-*^ arcsec^^ A~^. 

For flux calibration, we observed standard stars each night; the stars selected included 
EG 274, LTT 4364, LTT 6248, and LTT 7379, using equation 3 in §3.1 to obtain absolute 
fluxes. 

Most Fabry-Perot programs acquire observations of a diffuse monochromatic calibration 
source at intervals of 30 to 60 minutes during the night to monitor and correct the wavelength 
response. Our program yields high signal/noise ratio spectra of the atmospheric night sky 
lines; as described in §3 these can be used for wavelength calibration. Since they are observed 
simultaneously with the target, using them provides a more accurate calibration and also 
minimizes the need for additional calibration observations during night time observing. 

Daily calibrations were acquired. These included bias and dark frames, dome and twi- 
light flats with the etalon gap set to the spacings used during the night, and neon calibration 
cubes. We varied the tilt angle for the etalon from night to night, and obtained calibrations 
appropriate to each night. 
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Fig. 2. — (a) Spectrum for order 136. Two night-sky OH lines are labeled, along with the 
Galactic [N II] line. The spectrum is averaged over the entire field of view for which order 
136 is transmitted. The long wavelength cutoff occurs because the etalon gap was set so 
that the on-axis wavelength was ~6593 A. (b) Spectrum for order 135. (c) Blocking filter 
transmission percentage as a function of wavelength, normalized to the peak transmission. 
The transmission is averaged over the same region as for the spectrum shown in b). 
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3. REDUCTIONS 

Due to the narrow effective bandpass, low surface brightness and relatively small solid 
angle over which we expect a signal, we are severely shot-noise limited and in addition must 
deal with a number of subtle systematic effects which further degrade the signal. Therefore, 
in addition to the observing procedure described above, the reduction of the data must deal 
as well as possible with these difficulties and hence we describe in some detail various steps 
in the reduction of the data. We first give a simple overview of the basic reduction procedure 
and the approximations implicit in it for both the flux calibration observations and the target 
itself. 

3.1. Basic equations 

The count rate in ADU s^^ expected from a very small patch of sky (e.g., a pixel 
corresponding to solid angle AQ) may be expressed as 



{dc/dn)An = ge{\„,a.m A,jf{n,\,^a.m / Si\',n) R{\'-\^aA^)) d\' (2) 

In this expression S'(A, Vt) is the photon surface brightness per A as a function of wavelength 
and direction. At any given point on the detector, the Fabry-Perot will have its maximum 
response at some wavelength X^axi^)-, and at this wavelength and direction we denote the 
overall effective collecting area of the system as Aeff{^^\max{^)) (not including the atmo- 
spheric attenuation, e(A) and the CCD gain in ADU e^^, g, which are included separately 
in equation (2)). This effective area includes the spatial and wavelength dependence of the 
blocking filter and the pixel-to-pixel sensitivity of the CCD, determined via flat fleld images, 
as described below in §3.2. However, the etalon will have some response in this direction to 



-15- 

photons differing from Xmax{^) and we denote this response function relative to its maximum 
response by R{X — X^axi^))- In principle, as in any system with finite spatial and wave- 
length resolution, we would need to invert equation (2) to reconstruct the actual distribution 
of surface brightness, S{X, Vt). 

In practice, spectra for both the standard star and the target are constructed as follows. 
1) We define a fine grid in wavelength with a uniform spacing, AXgrid of 0.11 A. 2) For 
each pixel within the target mask (described below) we find the wavelength Xmax associated 
with the center of that pixel, calibrated using the series of neon lamp observations described 
below. 3) We associate all the counts in that pixel with the wavelength bin in this fine grid 
closest to Xrnax- After extraction on this 0.11 A interval grid, the spectra were then binned 
on a coarser 0.34 A grid, and all plots, fits and analyses of the spectra use the 0.34 A grid. 
We refer hereafter to each element of this coarser grid as a "wavelength bin" to distinguish 
it from the spatial pixels. 

Consider first the observations of an early type standard star for calibration of the 
throughput of the system. In this case, over the wavelength region of interest the fiux is 
constant with wavelength and the spatial structure is simply the seeing disk. In this case, 
the total count rate from the star will be 

Cstd = g e{X„iax{^)) Aeff{Xmax{^)) [dN/dX]std AXeb (3) 

In this expression, [dN/dX]std is the (known) monochromatic fiux (in photons cm^^ s~^ 
A~^), from which the effective area at this wavelength and this location on the detector 
may be calculated. The quantity AXeb is the effective bandwidth of the etalon, which is 
the response function R{X — Xmax) (equation 2) integrated over A. The response function 
of the Fabry-Perot is quite accurately represented by a Voigt profile, and we determine the 
effective bandwidth empirically using neon comparison lines. When the procedure described 
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above is applied to an unresolved continuum object, the only error is simply that associated 
with the aliasing of a detector which is not adequately sampled if the intra-pixel response is 
not uniform; however, this effect can be minimized by small dithering of the standard star 
observations. 

For the observation of the target, the opposite situation holds: the spatial scale is much 
larger than the pixel size, but the expected wavelength dependence in the neighborhood 
of the emission line is, in our case, comparable to the etalon resolution. In this case, the 
count rate from the source if it were finely sampled spatially can be shown to be just the 
convolution of the intrinsic emission line profile with the etalon response: 

/ + 00 
^(A' - Ao) i?(A' - A) dX' (4) 

-oo 

where B^ is the wavelength-integrated photon surface brightness, and G is the normalized 
intrinsic line profile. We have omitted the angular dependence for clarity. The instrumental 
profile i? is a Voigt function, and as discussed in §3.2 the convolution represented by the 
integral in equation (4) results in an observed profile which is also well described by a Voigt 
function. Then the photon surface brightness Bo is obtained from the count rate summed 
over the line profile F' = ^^ Fj AXgrid as 

F' 

^° = 4 A)— (^) 

g e Aeff AXeb 

The effective area A^/f = rj A, where A is the telescope area and rj is the efficiency 
of the entire optical path including telescope and instrument optics, blocking filter, etalon, 
and CCD. Using the standard star observations, we determined that rj = 0.22. Photon rates 
were converted to surface brightness, corrected for the small amount of Galactic extinction, 
which from Schlegel, Finkbeiner, & Davis (1998) amounts to about 0.05 mag in the R band. 
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In practice, as in conventional spectrographs, if the samphng due to the pixel size is not 
adequate, there may be aliasing of the spectrum. For the etalon order nearest the optical 
axis of the etalon, the sampling is adequate in spite of the 4x4 binning. However for the 
remaining orders this is no longer the case. If each wavelength bin were only sampled once, 
this could result in an emission line profile whose center and width were poorly defined. In 
fact, this situation is entirely mitigated because the format produced by the etalon causes 
each bin in the wavelength grid to lie at more-or-less random positions on the more than 
100 pixels contributing to each wavelength bin, just as "dithering" does in a conventional 
undersampled spectrograph. 

3.2. Details of the Reduction Procedure 

The first steps in the data reduction for all calibration and data CCD frames were bias 
and dark frame subtraction; bias and dark frames were constructed from the median value 
for each pixel derived from stacks of bias and dark frames in order to minimize the noise 
contribution from this step. Next, spatial variations in the instrument and CCD response 
were removed by dividing by a normalized flat-field frame obtained from a stack of exposures 
of the standard fiat-field lamp illuminating a screen on the side of the dome. We used a fiat- 
field frame obtained at the same etalon gap spacing / as used for the observation; we masked 
pixels for which the fiat-field response was less than 60% of the peak. 

The wavelength response given by equation (1) assumes ideal optics. In practice, there 
are deviations that result in slight shifts in the response; for optimum resolution we charac- 
terize these deviations as follows. We illuminate the dome screen with a neon lamp, stepping 
the etalon gap over a free spectral range in steps of half the wavelength resolution. Then 
at each pixel we determine the etalon gap that maximally transmits the neon calibration 
line. In this way, we characterize the deviations from the ideal wavelength relation given by 
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equation (1). 

As mentioned in §2.2, the standard method to monitor and cahbrate temporal variations 
is to take periodic exposures of a diffuse monochromatic cahbration lamp. Because we have 
bright night sky lines of known wavelength in each exposure, it is more accurate and efficient 
to use these. In principle, we could determine the centroid and radius of one of the night sky 
rings; however, in practice stars and cosmic rays make this difficult. Instead, we used the 
following method. The wavelength error caused by an incorrect choice of optical axis and 
etalon gap varies as a function of position in the field. The night sky lines present in the 
field (e.g.. Fig. 1) have known wavelengths. We compare the known wavelengths with those 
computed for a large range of assumed optical axes and etalon gap, and select the axis and 
gap which yield the smallest difference between the calculated and known wavelengths. 

We used the Tukey biweight statistic as described by Beers, Flynn, & Gebhardt (1990) to 
estimate the count rate per pixel for each wavelength bin based on all the pixels contributing 
to that wavelength bin. As described by Beers, Flynn, & Gebhardt (1990), the biweight 
estimator is more resistant against outliers and more robust than the mean, and more efficient 
than the median. It is highly effective at minimizing the effects of cosmic rays and stars. 

To check that variations across the frame do not cause offsets between the Ha part of 
the spectrum and the rest of the spectrum, we divided the frame into several regions and 
compared the spectra extracted from each region. Differences in these spectra were small, 
and are dominated by scattered light from stars; this effect is described next. 

Even though spectral features in stars in the field are likely to be weak, these stars can 
nonetheless create spectral artifacts. Recall that the count rate in each wavelength bin is 
determined from the biweight of the count rate for the CCD pixels corresponding to that 
wavelength bin. Although stars are effectively rejected by the biweight estimator, which 
excludes pixels with count rates that deviate by more than four standard deviations from 
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the central location, scattered light from brighter stars can cause difficulties. This is because 
portions of the extended skirts of the point spread function can be too weak to be rejected 
by the biweight estimator, yet extend over a sufficiently large number of pixels so as to bias 
the biweight. 

Accordingly, we adopted the following procedure. First, we excluded regions for which 
scattered stellar light was detectable by visual inspection of the frames. Since the error intro- 
duced by a star increases with the size of the region that escapes the biweight cutoff, and since 
the larger regions are of course associated with the brightest stars, this was straightforward, 
and in itself greatly reduced much of the problem associated with star light. This eliminated 
artifacts in individual spectra that were as large as 5 x 10^^*^ ergs cm^^ s^^ arcsec^^ A^^, 
although making a difference of less than 5 x 10^^^ ergs cm^^ s^^ arcsec^^ A~^ to the fi- 
nal averaged spectrum. As an added precaution, we identified the remaining stars in the 
field using the program SExtractor (Bertin & Arnouts 1996). Next, using a point spread 
profile determined on brighter stars, we subtracted the scattered light contribution from the 
remaining stars; the resulting changes to the spectra were minor. 

As described in §3.3, the incident ionizing flux is inferred from the Ha brightness toward 
regions with observed column densities N{HI) > lO^^cm^^ using the 75" H I maps provided 
in digital form by Chengalur, Giovanelli & Haynes (1995). Therefore, we excluded regions 
with lower H I column densities in constructing the spectrum for the wavelength interval 
expected for the redshifted Ha emission, 6588.9 A < A < 6591.9 A. Because pixels with 
wavelengths outside this range were not excluded, the number of pixels per wavelength bin 
over the wavelength interval defining the continuum was a factor of two to three higher than 
for the region of line emission. 

In Figure 3 we show the image resulting from directly adding the GHl and GH2 frames 
of one W-GH1-GH2-E set and subtracting the W and E frames. The OH night sky lines 
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largely cancel, while the [N II] line appears black, corresponding to brighter [N II] emission 
in the GHl and GH2 frames than the sky frames. Note the patchy spatial distribution of the 
[N II] emission and the presence of two well defined velocity components. Also note that the 
OH line residuals appear to have a slight offset in the subtracted image. This occurs because 
of drifts in the etalon. Because of such drifts, spectra are extracted not from W-GH1-GH2-E 
combined frames such as shown here, but from the individual frames because drifts in the 
wavelength scale can be corrected during the extraction process, as described above. The 
ring corresponding to the approximate location of the expected redshifted Ha signal with 
centroid at 6590.4 A is shown by the solid- line circle. 

Fig. 3. — Target minus sky frame produced from a single W-GH1-GH2-E set. Several OH 
and [N II] features are marked. Spatial variations in the Galactic [N II] features are evident, 
corresponding to differences in [N II] emission velocity and spatial distribution between the 
target and sky fields, along with OH night sky residuals. Also indicated are the N{H) = 10^^ 
cm~^ contour for the southwest HI 1225+01 cloud (GH SW) and for the bridge between the 
SW and NE components; the H I data are from Chengalur, Giovanelli & Haynes (1995). The 
cross marks the etalon optical axis. The solid circle shows the locus of 6590.4 A in order 
136, the expected wavelength for Ha, while the dashed arc shows the same locus for order 
135, where little or no Ha is expected. 

Target minus sky (i.e. On-Off) spectra produced for orders 136 and 135 from each 
W-GH1-GH2-E set of four spectra are shown in Figures 4 and 5; each On-Off spectrum 
corresponds to 30 minutes integration time on the target. The error bars shown on the bot- 
tom spectrum of Figures 4 and 5 are calculated from Poisson statistics; the errors are larger 
at the wavelengths of night sky lines due to the larger count rates, and in the wavelength 
range 6588.9 A < A < 6591.9 A due to the smaller number of unmasked pixels. Significant 
non-zero residuals are seen at the wavelengths of the Galactic [N II] line; this is due to the 
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difference in high Galactic latitude [N II] emission between the on and off positions. Also the 
atmospheric lines and the continuum level typically show non-zero residuals due to temporal 
and spatial variations in night sky brightness. 

Figure 6a is reproduced from Figure 2b to enable comparison with Figure 6b which 
shows the target minus sky spectra for orders 135 and 136 averaged over the 15 W-GHl- 
GH2-E spectra obtained during the run. In Figure 6c, we show a magnified view of the 
spectra, including error bars. The error bars in Figure 6c are estimated from the scatter at 
each wavelength between the 15 spectra that contribute to the average and have decreased 
by approximately a factor of four from the individual W-GH1-GH2-E spectra, as expected. 
Small negative residuals remain at the OH wavelengths, at a level of about 1% of the line 
strength, due to nonlinear temporal variations in the OH line strength. 

There is also a non-zero residual continuum level. We determined this residual con- 
tinuum level (assumed constant with wavelength over this narrow range) by averaging over 
6598.4 A < A < 6601.5 A and 6606.3 A < A < 6610.8 A for order 135. The residual contin- 
uum level is 1.4 x 10~^° ergs cm~^ s^^ arcsec^^ A^^, about 0.3% of the continuum brightness. 
As expected the percentage continuum residual is smaller than the OH residual since the 
OH lines are intrinsically more variable. 

There is a weak OH line at 6591.7 A evident in Figure 6a and which also has been 
identified by Osterbrock et al. from compilation of Keck spectra (Osterbrock, Fulbright, 
& Bida 1997). We fit each of the four OH lines in Figure 6a to determine the ratios of 
the OH line intensities. The 6591.7 A line is about 10% of the strength of the 6596 A 
line, similar to what can be derived from the Osterbrock, Fulbright, & Bida (1997) Keck 
data. Based on the residual negative OH emission of the bright OH line at 6596 A of 
about 2 X 10^^*^ ergs cm~^ s~^ arcsec"^ A~^, we expect a peak residual of approximately 
—2 X 10~^^ ergs cm~^ s^^ arcsec^^ A^^ for the 6591.7 line. We fit the residual fluxes in each 
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Fig. 4. — Target minus sky spectra for each of the 15 W-GH1-GH2-E sets obtained during 
the run are shown for order 136. The error bars shown on the bottom spectrum are calculated 
from Poisson statistics. 
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Fig. 5. — Target minus sky spectra for each of the 15 W-GH1-GH2-E sets obtained during 
the run are shown for order 135. The error bars shown on the bottom spectrum are calculated 
from Poisson statistics. 
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Fig. 6. — (a) Same spectrum as in Figure 2b showing the hues in the bandpass, (b) The 
average of all target minus sky spectra for order 136 (thick line) and 135 (thin line). Multiple 
components in the residual [N II] line are evident, due to differences in the spatial distribution 
of [N II] between the target and sky positions, (c) Expanded view of spectra shown in the 
panel above. Error bars are la derived from the scatter of the 15 W-GH1-GH2-E spectra, 
(d) The same spectra as in (c), except the continuum, OH, and [N II] residuals have been 
subtracted as described in the text. Here the error bars are those calculated from Poisson 
statistics. 
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of the OH lines shown in Figures 6b and 6c using a Voigt function; for the fit, the Gaussian 
and Lorentzian hne width parameters were held fixed using values determined from a fit to 
the neon calibration line. Also, the wavelength was fixed at the known value of the OH lines, 
and the continuum was fixed at zero since the fits were made to continuum subtracted data. 

The residual [N II] peaks in Figure 6b are at —50 km s^^ and +80 km s^^ respectively. 
The brightness of this [N II] feature is close to the detection threshold of the observations 
of known high velocity clouds made with the same instrument (Weiner, Vogel, & Williams 
2000), but with the sensitivity of the present observations it is a highly significant detection 
and clearly not an artifact. The brightness of the 80 km s~^ feature varies significantly on 
scales of roughly an arc minute, and arcs corresponding to this feature can even be seen in 
a few locations in Figure 1. 

In attempting to subtract off the [N II] from the spectrum it is necessary to model 
the feature with both components due to the broad Voigt wings in the etalon point spread 
function. Unfortunately, however, the strength and extent in velocity space of this [N II] 
makes the locating of the continuum around the expected position of the redshifted Ha 
more difficult, since the [N II] emission is clearly present to within ~ 140 km s^^ of the 
Ha and conceivably could extend to greater velocities at lower intensity. Moreover, the 
clear differences in structure of the [N II] between the on and off positions makes the on- 
off sky subtraction problematic over the spectral region around the [N II]. Note that the 
[N II] spectra even over the relatively small angles separating orders 136 and 135 are slightly 
different due to spatial variations in the [N II] . These can readily be seen by direct inspection 
of the difference images (e.g.. Figure 3). 

Although the presence of the [N II] adversely affects the measurement of the Ha signal 
which is our primary objective, it nonetheless is an interesting feature in itself. We are not 
aware of other observations in this general part of the sky that would have the sensitivity 
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to detect Ha or [N II] emission at the level we detect. In H I, the Leiden-Dwingeloo survey 
(Burton & Hartmann 1994) clearly detects emission only in the range —50 < Vlsr < +40 
km s~^. A similar result is found in the HIPASS survey (Barnes et al. 2001). Thus, we cannot 
at present associate the 80 km s~^ [N II] emission with known H I emission. The velocity of 
the 80 km s~^ feature is nearly large enough to qualify as a high velocity cloud (Wakker & 
van Woerden 1997); however, Wakker & van Woerden (1997) emphasize that the deviation 
velocity (the difference from the largest velocity permitted by Galactic rotation) is a better 
criterion, arguing for a deviation velocity greater than 50 km s~^ as the cutoff. Using the 
model presented by Wakker (1990), positive velocities due to rotation are not expected in 
this direction due to its high Galactic latitude {b = 64°); consequently the deviation velocity 
is 80 km s~^, suggesting that the detected emission may be produced by a high- velocity 
cloud. With an emission measure of 25 mR, the [N II] emission is comparable to the weakest 
Ha emission from known high velocity clouds observed by Weiner, Vogel, & Williams (2000); 
however, [N II] emission was not detected by Weiner, Vogel, & Williams (2000) toward the 
faintest high velocity clouds, with upper limits somewhat below our detection. In other words 
the [N II] emission we detect is somewhat brighter than in some known H I high velocity 
clouds, in spite of the fact that no H I emission is seen in this direction at this velocity. 
The detected [N II] emission is clearly interesting and merits further study to determine its 
origin. 

The resulting spectra with the continuum, OH line, and [N II] line residuals subtracted 
are shown in Figure 6d. In this figure, the error bars have been calculated from Poisson 
statistics. Note that the order 135 and 136 spectra are generally the same within the sta- 
tistical errors, except for part of the [N II] line; this is due to spatial structure in the [N II] 
distribution and the fact that the difference in [N II] distribution between the target and sky 
positions must be slightly more complex than can be modeled with just two components. 
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Next we fit a Voigt profile to this spectrum over the wavelength interval 6588 A < A < 
6593.5 A (corresponding to about four spectral resolution elements). The continuum is fixed 
to be zero since the continuum has been subtracted, as discussed above. The instrumental 
profile is taken from the neon calibration. 

The assumed intrinsic Ha line profile was obtained as follows. Using the 75" resolution 
VLA H I observations of Chengalur, Giovanelli & Haynes (1995) (which are given in the 
heliocentric frame), we determined the H I velocities of the gas at each pixel on the Fabry- 
Perot ring at the approximate redshifted wavelength of Ha and then shifted these velocities 
by +5.8 km s~^, which puts them in the barycentric frame which we have used throughout 
this work. Weighting each pixel equally as a first order approximation since the Ha emission 
is only weakly dependent on H I column density for the observed column density range, we 
constructed the spectrum shown in Figure 7a. H I emission occurs mostly at two velocities 
because HI 1225+01 has a rotation curve characteristic of a galactic disk, and the 6590.4 
A ring has relatively little overlap with the central part or minor axis (e.g.. Fig. 3). This 
profile was then convolved with a Gaussian with a dispersion of 6 km s^^ (corresponding to 
a typical dwarf dispersion and slightly less than the velocity dispersion of 8 km s^^ of the NE 
component) and also with a Gaussian whose width corresponds to the thermal velocity of 
the ionized hydrogen at a temperture of 10000 K. This produces the profile shown in Figure 
7b. Finally, this profile was convolved with the instrumental profile to produce the profile 
shown in Figure 7c. This corresponds to the profile of the Ha line we expect to observe from 
HI 1225+01. This is well fit by a Voigt profile centered at 6590.4 A and FWHM 79 km s"^ 

A Voigt profile with this shape is then used as the line profile in a fit to the observed 
spectrum, in which the only free parameter is the Voigt line strength. The best fit line 
brightness is 4.2 mR. The spectrum and the 4.2 mR best fit for the OH, [Nil], and continuum 
subtracted spectrum for order 136 are shown in Figure 8b. To estimate the uncertainty in 
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Fig. 7. — (a) Wavelength distribution of redshifted Ha centroids expected based on 75" VLA 
maps of H I obtained by Chengalur, Giovanelli & Haynes (1995) for Ha emission regions 
observed with the Fabry- Perot, (b) Expected intrinsic Ha profile obtained by convolving the 
H I profile in (a) with Gaussians for H atoms at 10000 K and the H I velocity dispersion, (c) 
Expected observed Ha profile, obtained by convolving the intrinsic Ha profile in (b) with 
the Fabry-Perot instrumental profile. 
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the line strength, we used a Monte Carlo approach (Weiner 2001), in which we constructed 
spectra using the errors estimated from the comparison of the individual W-GH1-GH2-E 
spectra (i.e. the errors shown in Figure 6c). The resulting la error bars are ±1.6 mR, and 
the 2a upper limit is 7.4 mR. 

Although we have attempted to accurately apply flatfield corrections to the data, differ- 
ences in the spectral and spatial properties between the flatfield lamp and the sky data, as 
well as differences between the resolution in orders 135 and 136 mean that direct subtraction 
of the "on-only" data in order 135 from 136 to obtain the net signal in the region of the 
expected Ha signature may well mask the signal we are attempting to detect. For this reason 
we have employed the on-off chopping described previously, which should accurately remove 
the effects of these differences in the instrumental response. 

However, although the sequence of on-off observations removes these differences as well 
as temporal variations in the night sky background which are linear, the characteristic time 
for the full cycle for one off-on-on-off sequence is over an hour, and variations over a shorter 
time may well occur. Therefore the order 135 spectrum of Figure 6d was subtracted from 
order 136 in the same figure to produce Figure 8c. (In order to reduce the noise in this 
subtraction the data in order 135 were first smoothed using a 5-point running mean.) 

Disappointingly, the signal present in Figure 8b appears to have largely disappeared. 
It is known that other weak night sky emission features are present. In particular, weak 
emission from O2 has been found in the night sky spectrum (Slanger et al. 2000) and one 
of these emission features falls at the expected wavelength of Ha (T. Slanger 2001, private 
communication). No attempt has been made to model any residual emission from this 
source. However, the impression is that the main effect of the subtraction of 135 from 136 is 
simply to lower the level of the continuum so that it falls below the nominal zero level. This 
impression is confirmed by again fitting a Voigt profile to the same data as shown in Figure 
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8c, but now allowing the continuum level to be a second free parameter. The resulting best 
fit, shown in Figure 8d, yields 3.9 mR, very similar to that found above. Because of the 
very limited continuum region in order 136, which is affected nearly up to the Ha signal at 
shorter wavelengths by the [Nil], as described above, and only extends to ~ 6593 A which 
is the termination of order 136 as explained in §2.2, it is not clear whether such an offset 
really exists. 

It is also not clear what could give rise to such a continuum offset. Spatial gradients 
in the background continuum arising from the zodiacal light and airglow continuum^ would 
have to be present over an extent of only ~8'. It should be borne in mind however that 
the possible offset is only about 0.1% of the total continuum level. The data provide no 
evidence for such gradients; first, the CCD dewar was rotated 180° after the second night, 
yet comparison of spectra obtained from the opposite rotations does not show mirrored 
spectra artifacts; second, comparison of the order 135 and 134 spectra in Figure 9d reveals 
no significant difference, even though these are extracted from regions offset by ~ 5', nearly 
as large as the 7-10' offset between orders 136 and 135. 

A more likely possible source of artifacts is scattered light from stars that has not entirely 
been removed or masked. This scattered light would be more likely to create spectral artifacts 
in the central order. Partly this is because out of focus reflections between the etalon and 
blocking filter are not directed out of the beam for the central order due to the etalon tilt 
angle selected. Also, due to the quadratic nature of the wavelength dispersion scattered star 
light affects relatively narrow spectral regions in the central order compared to the outer 
orders. This may explain why the order 135 and 134 spectra (e.g.. Figs. 8a) are quite 
similar, and more featureless than the order 136 spectrum. 



^This airglow is generally ascribed to NO2, although there is some disagreement on this point. 



-31- 

Although there is some weak indication of the detection of a hne from Figures 8b and 
8d, with strengths of ~ 4 mR, because of this uncertainty we cannot claim a positive direct 
detection of the ionizing background radiation, but the formal 2a upper limit^ to any signal 
from the fits in either 8b or 8d is about 8 mR, and we shall use this limit in our subsequent 
discussion. The actual formal upper confidence limits to the fits in Figures 8b, c, and d are 
shown in Figures 9b, 9c, and 9d, where Figure 8a is again reproduced for reference as Figure 
9a. 

3.3. Conversion of the Ha Surface Brightness to Ionizing Flux 

Following VWRH and Stocke et al. (1991), (see also SRGPF) the one-sided incident 
ionizing flux is defined as 



oo 



$ = TT / ^ rfz/ (6) 

where J^ is the mean intensity. The factor tt assumes that the background intensity is 
isotropic. From our measurement of the Ha surface brightness, I Ha-, we may calculate $ 
from: 



-' Ha -L / /ipj,QJ 



hi' JajHa V Atot 



which can be written as 



^The Monte Carlo calculations actually lead to a direct estimate of an upper confidence limit of 97.5%, 
which for a Gaussian distribution would imply that 95% of the Monte Carlo trials lie between -2cr and +2cr 
of the best fit value, but we will simply refer to this as the "2cr" upper limit. 
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Fig. 8. — (a) Target minus sky spectra compared for orders 136 (thick solid line), 135 (thin 
solid line), and 134 (dashed line). In this panel, unlike those below, OH and [N II] residuals 
have not been subtracted, (b) Order 136 spectrum with continuum, OH and [N II] lines 
subtracted, shown as a thin line; the error bars are estimated from the data. The thick line 
is a fit of a Voigt profile in which the only free parameter is the line strength. The best 
fit emission line has strength 4.2 ±1.6 mR. (c) As in (b), except a smoothed order 135 
spectrum has been subtracted. The Voigt fit has strength 0.8 ±1.6 mR. (d) As in (c) except 
that the fit has both line strength and continuum level as free parameters. The best fit is 
3.9 mR ± 2.2 mR. The best-fit continuum is shown as a horizontal dashed line. 
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Fig. 9. — (a) Target minus sky spectra compared for orders 136 (thick solid line), 135 (thin 
solid line), and 134 (dashed line). In this panel, unlike those below, OH and [N II] residuals 
have not been subtracted, (b) Order 136 spectrum with continuum, OH and [N II] residuals 
subtracted, shown as a thin line; the error bars are estimated from the data. The thick line 
shows a 7.4 mR emission line, corresponding to the 2o" upper limit, (c) Order 135 spectrum, 
with 7.4 mR emission line overlaid, (d) Same as (c), except the error bars are Poisson. 
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$ = (1.0 X 10=^ cm-Vi) f^) —^ (^\ (8) 

^ ' KmRJ fafHa yAtotJ ^ ' 

In equations (7) and (8), fa is the fraction of the one-sided incident flux which results in 
recombinations to excited states, fna is the fraction of those recombinations which produce 
Ha photons, Atot is the total surface area exposed to the incident radiation and Ap^-oj is the 
observed projected surface area. This expression is quite general, except that for complex 
geometries the appropriate averages over the cloud surface of the quantities fa and fua need 
to be taken, and the lua surface brightness is that averaged over the projected area. 

3.3.1. Effective Ha Recombination Rate 

The Ha spectrum has been extracted from regions for which Nhi > 10^^ cm~^. However, 
as described in section 3.3.2, recent work by Chengalur, Giovanelli & Haynes (1995) suggests 
that HI 1225+01 SW is a disk which is viewed far from the normal, with an aspect ratio 
of 4:1, so that the ratio of the total area to the projected area is ~ 8. If for simplicity we 
assume plain slab geometry, then the column density normal to the disk corresponding to 
the lO^^cm"^ column density is 2.5 x lO^^cm"^, corresponding to an optical depth at the 
Lyman limit of 15.7. Adopting the spectral index of 1.8 suggested by SRGPF we find that 
at this optical depth about 87% of the normally incident photon fiux is absorbed within 
the disk. However, some of the resulting recombinations occur to the ground state, and a 
portion of this diffuse Lyman continuum radiation can escape from the slab after one (or 
more) such scatterings, and thus does not contribute to any Balmer radiation. The fraction 
of the recombinations which go to the ground state depends upon the electron temperature, 
as does the energy distribution of this recombination radiation, and as discussed below, we 
adopt an electron temperature of 10000 K. Since most of the region around the Fabry-Perot 
6590.4 A ring consists of H I column densities in excess of 10^^ cm"^, we calculate the average 
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value of fa by creating a histogram of column densities which contribute to the ring (after 
reducing the observed values by the geometric factor of four) and for each of the values of 
the column density in the histogram compute the value of fa, weighting each of these values 
by the number of pixels contributing to that bin. Detailed calculations then yield a fraction 
of the normally incident photon flux which ultimately results in recombinations to excited 
levels of 0.88, which is the quantity fa in equations (7) and (8). (Had we used a spectral 
index of 1.0, the number would be reduced to 0.75. We also note that even for a normal 
column density 10 times lower than our minimum value the fraction of the continuum flux 
which results in excited state recombinations is still slightly greater than half of the 0.88 
value, namely 0.42.) 

We also require the effective recombination rate for producing Ha photons, or in the 
formulation used in equations (7) and (8), the fraction fna of excited state recombinations 
which ultimately produce an Ha photon. The value of fna depends, as does fa, upon the 
electron temperature, and upon whether Lyman line photons readily escape following the 
recombinations ( "case A" ) or scatter internally until they are ultimately converted to Balmer 
and higher series photons ("case B"). Since nearly all of the recombinations occur at optical 
depths in the Lyman lines which are very large, case B is appropriate. 

In VWRH we adopted a temperature of 20000 K on the assumption that a hard spectrum 
and low metallicity would lead to an electron temperature much higher than the value of 
~ 10000 K usually assumed. However this assumes a balance between atomic heating and 
cooling mechanisms. In regions of the low density intergalactic medium producing the low 
column density Lya forest lines, Dave & Tripp (2001) deduce much lower temperatures 
(~ 5000 K) due to the fact that Hubble expansion dominates over atomic cooling processes. 
However, these considerations may not be applicable to the material comprising HI 1225+01 
where the expansion has evidently been halted, and where the particle density is probably 
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much higher than that in which the low column density Lya lines arise, in which case thermal 
equilibrium is a much better approximation. Thermal equilibrium calculations for a range 
of metallicities and particle densities for gas exposed to metagalactic radiation and partially 
shielded by neutral H I yield temperatures in the range 8000-10000 K (Wolfire 2001 private 
communication). However his calculations do not extend below column densities of 10^^ cm^^ 
so that on the outer skin of the cloud the temperatures might be slightly higher. We adopt 
an electron temperature of 10000 K. Then, for case B and an electron temperature of 10000 
K from data in Osterbrock (1989) and Pengelly & Seaton (1964) we find fna = 0.45. 

3.3.2. Cloud Geometry 

As noted previously, Chengalur, Giovanelli & Haynes (1995) have made a detailed study 
of the H I in HI 1225+01 using the VIA with 40" resolution. They deduce that the SW 
component is a thin rotating disk viewed nearly edge on, with an axial ratio of about 4:1. 
This is based upon the shape of the contours of the higher column density regions, but also 
on the shape of the velocity displacement as a function of position as well as the fact that 
such a steep inclination would explain why the peak column density is higher than in the 
NE component. They also point out that if the aspect ratio is > 3 then the true (normal) 
surface density in the SW component is not high enough to trigger star formation (none 
is observed), contrary to the situation in the NE component. For this proposed geometry, 

^tot/^proj = 8. 

However, the dynamical mass for HI 1225+01 SW calculated by Chengalur et al (1995) 
within a radius of 3' is 7 x IO^Mq, slightly smaller than the total HI mass of HI 1225+01 
SW. An HI mass comparable to the dynamical mass is of course contrary to what is now 
known about the ratio of total to baryonic mass. This apparent discrepancy can be reduced 
in several ways. First, the HI mass refers to the total HI mass of HI 1225+01 SW while the 



-37- 

dynamical mass is only calculated for the inner 3'. Second, the distance to HI 1225+01 may 
be as much as a factor of three closer than 20 Mpc, the distance assumed by Chengalur et 
al (Salzer et al. 1991; Impey et al. 1990). This would reduce the HI mass more than the 
dynamical mass. Nonetheless, the rotational velocity of HI 1225+01 SW, 13 km s"^, with 
just a small correction for inclination if the Chengalur et al geometry is correct, is extremely 
low for the HI mass, even at the near distance — see the Tully-Fisher type relation between 
rotational speed and HI mass derived by McGaugh et al. (2000). Indeed, a reasonable 
interpretation of the observed morphology of HI 1225+01 SW is that it has two spiral arms 
elongated along the line towards HI 1225+01 NE, with a more face-on inclination. 

Accordingly, we must conclude that the geometry is uncertain. The extremes are repre- 
sented by the 4/1 inclined disk, with Atot/Aproj = 8, and a face-on disk, with Atot/^proj = 2. 
For comparison, a sphere gives Atot/Aproj = 4. An infinitely long cylinder with the long axis 
in the plane of the sky, gives Afot/Aproj = tt; this factor is larger if the cylinder is inclined to 
the line of sight. An ellipsoid rotated about its major axis and with a major to minor axis 
ratio of 4 (perhaps a reasonable model for a tidal feature) gives Atot/Aproj ~ 3.22. Again, 
this factor would be slightly larger if the major axis is inclined somewhat to the line of sight. 
In the discussion in §4 We adopt Atot/Aproj = 4 as a fiducial factor, recognizing that there 
is an uncertainty of about two in either direction. We will refer to the Atot/Aproj values of 
2, 4, and 8 as "low", "fiducial", and "high" in the following discussion. 
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4. RESULTS AND DISCUSSION 

4.1. Inferred limits on the Radiation Field and Comparison with Previous 

Results 

Using the values for the geometrical factors just discussed, the 8 mR upper limit for 
the Ha flux, and values of fa and f^a from equation (8) we obtain values for the one sided 
ionizing flux of $ = 10.1, 5.04, and2.52 x 10^ photons cm^^s^^ for the low, flducial and 
high values of the geometrical correction factor Atot/Aproj- The flrst two of these are within 
the range inferred by Maloney (1993), while the third is about a factor of two lower, though 
it should be remembered that our surface brightness value is an upper limit. Expressing 
the radiation fleld as J^ = J^J^^{v/vQ)~°', we obtain J,^,, = a7r^^/i$, where h is Planck's 
constant and a is the spectral index. Assuming the SRGPF value for the spectral index 
of a = 1.8, the corresponding upper limits for the mean intensity at the Lyman limit are 
Jj,(, = 0.96, 1.92, and3.84 x 10^^^ ergs cm^^ s^^ sr^^. The local unshielded metagalactic 
photoionization rate is calculated from this as F = At^: J^J^^av^^h^^ {a + 3)^^, which yields 
log(r) = —13.62, —13.32, and — 13.02, where F has the dimensions s^^. Had we adopted 
a different spectral index (e.g., 1.0) the corresponding numbers for our flducial value of 
Atot/Aproj = 4 of $, J^„ and log(F) would be 5.92 x 10^ 1.24 x 10"^^ and -13.43, respectively. 
If our possible detection at 4 mR should ultimately be conflrmed, the actual values would 
of course all be lower by a factor of two. 

A review of recent attempts to detect, or set limits on, the intensity of the ionizing 
background may be found in SRGPF. Some of these are summarized in their Table 1. It 
should be noted however that in the paper by Donahue, Aldering & Stocke (1995), their 
deflnition of $ is in terms of the ionizing speciflc intensity, not the flux. They also assume 
a value for fa of 1.0, so their numbers in their Table 2 cannot be compared directly with 
ours. Their upper limit on the Ha surface brightness from HI 1225+01 SW is about 65 mR. 
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The Ha surface brightness limit of VWRH (20mR) is the lowest surface brightness limit we 
are aware of prior to the current work, which is about 2.5 times higher than our new upper 
limit. Our use of different assumptions about the values of fa and f^a and the range of the 
projection factors considered here (VWRH assumed a face-on slab) results in our current 
value for $ which ranges between about 2.5 to 10 times lower than that quoted in VWRH. 

4.2. Comparison with Estimates of AGN and Hot Star Contributions to the 

Background 

A recent careful analysis of the contributions expected from both AGNs and hot stars 
to the ionizing background has also been carried out in SRGPF as well as estimation of 
the uncertainties in these estimates. In principle, these contributions are directly observable 
independent of cosmological parameters, but of course local H I prevents radiation near the 
Lyman limit from penetrating into the disk of the galaxy along essentially all lines of sight 
from the sun. In practice, according to SRGPF, the two largest uncertainties in the estimate 
of the contributions from the AGNs are the UV luminosity functions for the AGNs and the 
assessment of the amount of attenuation from the cumulative absorption of intervening H I 
between ourselves and the AGNs. 

As noted above, SRGPF derive a best estimate for the UV spectral index of 1.8 which 
we have adopted. For the value of the specific intensity at the Lyman limit they quote 

J.o,agn = l-^t^ol >< lO^'^^ crgs cm-^s^^ Hz^^sr-i (9) 

For the contributions from hot young stars embedded in galaxies the luminosity function 
of ionizing radiation is also somewhat uncertain, but the dominant uncertainty in this case is 
the fraction f^sc of ionizing photons which escape the galaxy. SRGPF thus characterize this 
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fraction as a free parameter but favor a value of order 5%, and for this contribution quote 

Ju,,stars = l.llaT >< lO^'^^ ergs cm-\s~^ Hz~\sr~^ (10) 

the sum of which is to be compared to our 2a upper hmits of J^,„ which range from 
0.96 to 3.84 X 10"^^ ergs cm^^s"^ Hz^^sr-^ 

If we examine the various combinations of the "best" , "upper" and "lower" la estimates 
in SRGPF we find that our largest value (for the face-on projection) is compatible for all 
but the largest of these combinations. For our fiducial value of the projection factor, our 
limit is compatible with, but slightly lower than, the sum of the two "best" contributions. 

4.3. Relevance to the Lya Simulations 

The range of our upper limits on the inferred z ~ photoionization rate of log F = 
— 13.6 to — 13.0 may be compared with that inferred from very recent numerical simulations 
of the Lya forest by Dave & Tripp (2001). 

The recent advances in both computing power and numerical simulation techniques have 
produced results for the uncondensed gas giving rise to the Lya forest which reproduce with 
impressive faithfulness the observed properties of the Lya absorption lines over the entire 
range of observed redshifts from ~ to ~ 5. Assuming a specific cosmological model, it is 
possible to infer, at any epoch, the value of the ionizing flux which can be compared with 
empirical determinations such as the one we have presented above. 

The basic argument is as follows: A complete specification of the cosmological param- 
eters (e.g., Hq, Qbary, ^A, the H/He ratio, and, e.g. an LCDM initial perturbation model) 
completely specifies the total baryon density at any redshift, and, in principle the course of 
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the entire process of gas condensation. The simulations thus yield estimates of the actual 
density of hydrogen nuclei at any point in the spatial grid and at any redshift, along with 
the velocity field of that material, and thus the local column density of hydrogen nuclei. 
Comparing these simulations with the observed values of the frequency of the H 1 absorption 
lines one can then infer the required fraction of hydrogen nuclei which are neutral and thus 
the H I photoionization rate, F. Dave & Tripp (2001) deduce a value of logF = —13.3 ± 0.7 
at a mean redshift of 0.17, and a rough extrapolation of this estimate to zero redshift yields 
a value of about —13.5, which is compatible with our range of upper limits. Of course there 
are considerable uncertainties in the simulations themselves, in addition to the basic cos- 
mological parameters, but these are most severe for the highest density regions where the 
relatively rare higher column density H I systems are formed. With further improvements 
in the sensitivity of the Ha measurements, improved data on the density of the low redshift 
Lya forest line density, and refinements in the simulations themselves, it will be possible to 
impose a constraint on the combination of the cosmological parameters which are involved 
in the Dave & Tripp (2001) estimate, especially those that yield the present epoch baryon 
density. 

4.4. Caveats 

As our detailed discussion of the reduction procedures and our final spectra make clear, 
there are certainly uncertainties in our fitted value for the Ha surface brightness and its 
upper limit. However, the uncertainty in the geometry of HI 1225+01 SW also introduces 
an additional major uncertainty. 

There is also the possibility that our (implicit) assumption that the H I gas has a pro- 
jected filling factor near unity is wrong. Suppose in fact the H I gas in the cloud has already 
undergone significant fragmentation on scales much smaller than the 21 cm resolution, with 
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each "cloudlet" having H I column densities substantially higher than that inferred under the 
assumption that no such small scale structure exists. This implies a corresponding reduction 
in the actual solid angle subtended by the high column density fragmented gas. In such 
circumstances the surface brightness of Ha from each little cloudlet would be only slightly 
increased from that expected in the absence of fragmentation, however the total Ha flux 
would be reduced by the projected area filling factor, and the true value of the photoion- 
ization value would be higher than that we have inferred. We consider this scenario rather 
unlikely since one would expect star formation to be underway in at least some portions of 
the structure if such gas fragmentation had already occurred to such an extent. It is also 
possible that the tidal interaction between the NE and SW components has produced fila- 
mentary or clumpy HI structure as sometimes seen in higher linear resolution observations 
of tidally interacting galaxies. 

We have also considered the possibility of significant amounts of internal extinction of 
Ha in the SW component by dust. Even for a normal dust/gas ratio, for N{HI) = 10^^ 
cm^^ the extinction is negligible; however, for the peak observed column density of 1.1 x 10^^ 
cm~^ in the 40" resolution map of Chengalur, Giovanelli & Haynes (1995), the emission in 
that direction could be attenuated to 75% of the extinction-free value. However, most of the 
6590.4 A ring lies toward regions of the target where the column density is much lower (see 
Figure 3), with a median column density of 4 x 10^^ cm^^. In any case, a normal dust /gas 
ratio in a cloud in which no star formation seems to have occurred seems unlikely. 

It would be useful to reobserve the SW component of HI 1225+01 with improved sensi- 
tivity in order to see if Ha at the level of ~4 mR is really present. However, the ambiguity 
in the geometry of HI 1225+01 SW remains a difficuly. It would be very interesting for 
numerical simulations of interacting galaxies to be run to see if the morphology and velocity 
fields of the entire HI 1225+01 complex can be understood, and this may shed some light on 
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the correct geometry. Beyond the issue of the geometry however, investigation of the initial 
conditions and subsequent evolution of the rare large isolated H I clouds such as HI 1225+01 
in which little or no star formation has taken place would seem a worthwhile program. In 
any event, it will be important to observe with increased sensitivity several other suitable 
targets, since although the SW HI 1225+01 component has some features which make it an 
attractive target, the [N II] Galactic emission and the uncertain geometry make a definitive 
determination of the background flux problematic. 

The expert assistance of Bill Kunkel and Oscar Duhalde in resolving technical problems 
at the telescope is gratefully acknowledged. RJW thanks Romeel Dave for helpful discus- 
sions on the Lja forest simulations and Don Osterbrock and Tom Slanger for enlightening 
discussions and data concerning the night sky spectrum. SNV and SV acknowledge support 
from NSF AST-9529167. Finally, we acknowledge the helpful comments of the referee, P. 
Maloney, in particular for pointing out the dynamical objection to the Chengalur, Giovanelli 
& Haynes (1995) geometrical model. 
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